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Abstract. We investigate the evolution of a central spin coupled to a spin bath 
without internal dynamics. We compare the cases where the bath couples to one or 
two components of the spin. It is found that the central spin dynamics is enhanced in 
the latter case, which may be interpreted as a frustration of dissipation. However, the 
quantum purity of the spin decays fast in both scenarios. We conclude that symmetric 
coupling of the bath to two orthogonal components of the spin inhibits dissipation but 
not decoherence. 
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1. Introduction 

The study of the role of a dissipative environment is of central importance to the field 
of quantum computation and for the fundamental understanding of the transition from 
quantum to classical behavior. In that context, the dissipative two-level system (TLS) 
is a well-studied paradigm [U [2]. The generic TLS, which in particular can be a spin | 
particle, experiences dissipation due to its coupling to a bath of harmonic oscillators. In 
the resulting spin-bath problem, an external magnetic field interacts with one component 
of the spin operator (e.g. S^) while a second component (e.g. S^) couples to an oscillator 
bath. Depending on the relative strength of the two interactions, one can switch over 
between underdamped and overdamped behaviour, as may become manifest through 
quantities such as the average spin energy or various correlation functions. Recently 
Castro Neto et al. [3, 4] have shown that if two components (5*"^ and S*^) of the 
effective spin are coupled to two different baths, then the competing effects of those 
baths can reduce the effect of dissipation. In particular, they have found that, for a given 
coupling strength, symmetric coupling to the two spin components is less decoherent 
than coupling to a single component. Moreover, in the case of symmetric coupling to 
and S*^, coherent behaviour is preserved for arbitrarily strong coupling. These two 
properties are remarkable because one would naively expect that coupling to a higher 
number of bath oscillators would increase the effect of dissipation. Rather, they have 
shown that the competition between the baths contributes to protect the TLS energy 
gap. The reduction of dissipation stems from the non-commutative character of the 
spin operators coupled to the different baths. Logically, if the two baths interact with 
the same spin component, then no dissipation reduction is observed. In Refs. [3],|1], it 
is argued that this feature arises from the lack of a preferred basis to which the TLS 
may relax at long times. This new phenomenon has been coined quantum frustration of 
decoherence, since it is interpreted as the frustrated attempt of the two environments 
to "measure" simultaneously two non-commuting observables. The study of Refs. [3], H] 
has been restricted to equilibrium properties such as the transverse susceptibility, which 
was evaluated using the numerical renormalization group. 

Features of quantum frustration were also reported in Ref. [5], where it was noted 
that the phase-number variable of a superconducting Josephson junction is coupled 
simultaneously to two different dissipative environments: through the phase to the 
quasiparticle field and through the Cooper pair number to the quantum electromagnetic 
field. The result is that the uncertainty in the macroscopic phase has contributions from 
both baths which tend to cancel each other. They never cancel completely because, in 
that particular case, the sources of dissipation differ widely both in nature and strength. 
Frustration of decoherence in Josephson networks has also been investigated [6] 

In Refs. [3, E] the dynamics of an oscillator coupled through its position and 
momentum to two different oscillator baths was studied. The problem can be shown 
to be equivalent to that of a large (quasiclassical) spin impurity in a ferromagnetic 
environment. It was noted that the dissipative oscillator may be driven from overdamped 
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to underdamped behaviour by the symmetrical addition of a second bath. This 
surprising effect was found only for the case where the momentum operator [p) and 
the position operator (g) are coupled to two different baths. Like in Refs. P H] , these 
effects were investigated through equilibrium properties such as the position-position 
response function. However, some dynamical aspects were also analyzed, noting that 
the purity of the quantum oscillator decays faster in the presence of two baths than 
in the presence of a single bath. This last result is important because it reveals that 
the competition between two baths coupling to non-commuting observables is not a 
universal panacea to suppress decoherence. Our present study is motivated by the need 
for a more detailed understanding of the equilibrium and dynamical properties of a 
quantum system in the presence of competing environments. 

In the studies of quantum frustration made earlier P HI [5], El [8] , the bath has been 
modeled by a set of non-interacting oscillators. However, if the interpretation is correct 
that the essence of quantum frustration stems from the canonically conjugate character 
of the two observables which couple to separate baths, then one should expect a similar 
behaviour to appear when the dissipative environment is formed by a bath of spins 
acting on a central spin impurity. We note that, because of the vector nature of the 
spin, the spin bath can by itself be viewed as formed by several baths. Thus a single 
spin bath may exhibit features of quantum frustration. Spin baths have been studied 
P, [10, [m [El [131 [m [El [E] as an alternative to the conventional oscillator models of 
quantum dissipation [H [2]. They are known to give rise to non-Markovian evolution, 
with the system evolution showing a strong dependence on the polarization of the initial 
state [15j . 

In this work we shall evaluate the dynamical properties of a TLS coupled to a bath 
of spins. Similar to the work of Novais et. al. [31 H] , we shall consider the situation where 
the components of the spin are coupled to two different baths and study the effects of 
frustration arising due to the non-commuting nature of the spin components. The work 
of Refs. [31 [Ij was based on the method of the numerical renormalization group. By 
assuming that the bath has no internal dynamics of its own, we are able to perform 
an analytical study. Comparing the case of coupling to a single component to that of 
symmetric coupling to both components, we find that the spectral function behaves 
similarly to the study of Refs. [31 H]. Namely, the spectral function develops a peak in 
the symmetric case which is absent in the case of single-component coupling. In Refs. 
[31 [1] , this was interpreted as the preservation of decoherence arising from the frustrated 
attempt of the two environments to measure two non-commuting observables. However, 
we find that the emergence of the peak as the second bath intervenes is compatible 
with a fast decay of the quantum purity of the central spin. This means that the 
frustration induced by the two competing environments has more to do with dissipation 
than with decoherence. Energy relaxation is indeed inhibited by the presence of a second 
bath coupled to the other spin component, while the quantum purity decays fast in the 
presence of a second bath, with only a minor form of frustration revealed by a short-lived 
revival which will be discussed. 
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Section 2 is devoted to the presentation of the model of a central spin coupled to 
a spin bath. In Sec. 3 we study the time evolution analytically, deriving expressions 
for the expectation values and the quantum purity of the central spin. Sections 4 and 
5 focus on the density of states and the response function of the central spin. Section 6 
deals with the tailoring of the spin bath properties which mimics the behavior of a spin 
coupled to a conventional bath of harmonic oscillators. Finally, the main conclusions of 
this work are summarized and discussed in Sec. 7. 

2. Central spin model 

We consider the dynamics of a two-level-system which is linearly coupled through two 
non-commuting observables to two independent environments of two-level-systems. We 
will refer indistinctively to both the central impurity and the constituents of the bath as 
particles of spin ^ or two-level-systems. If the spins of the environment carry their own 
dynamics, in general there are no conserved quantities other than energy and the system 
can not be treated analytically without approximations. However, in most of the solid- 
state spin systems where the spin bath interaction is a dominant source of mechanism for 
the dissipation of a TLS, the internal bath dynamics is generally very slow (for example 
in quantum dot systems, where the bath spins are nuclear spin half particles and the 
TLS is the electronic spin, [H]). We therefore assume that both environments carry no 
dynamics of its own, i. e. that their Hamiltonians are zero. Thus the bath dynamics is 
exclusively due to its interaction with the central spin [16]. The total Hamiltonian of 
system and bath is then given by 

H = Hs + HsB (1) 

where 

Hs =cooS', (2) 

N N 
k=l 1=1 

where S"*, i = x,y,z are the components of the spin operator of the central spin and 
and J| are spin operators of the bath spins. We assume homogeneous interaction 
between the central spin and the baths. Moreover we assume that the number of spins 
of each environment to be the same. The strength of the coupling to each environment 
is thus described by one parameter Qi [i = 1,2) only. For a single bath-environment 
the case of non-homogeneous coupling {gi, g2 dependent on index k, I, respectively) was 
solved explicitly by in Ref. [10]. However, it was shown in Refs. [lOl [16] that the 
only effect of inhomogeneous interaction is that of destroying certain revival effects. 
All other features can be captured within the homogeneous interaction between system 
and bath. The main advantage of the homogeneous interaction approximation is that 
exact, closed- form expressions for the expectation values and correlation functions can 
be obtained. 
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The Hamiltonian ([T]) is similar to that employed in Refs. [21 S] in that two different 
environments couple to the two perpendicular components of the central spin. The 
main difference with the model of Refs. [31 H] is the non-dynamic character of the 
bath, which we consider, which contrasts with the oscillator bath there considered. We 
shall examine whether our simpler model ([T]) can yield frustration effects similar to 
those obtained from the more complex model of Refs. [31 H], which was solved with the 
numerical renormalization group method. 

In the particular case where the spin bath couples to only one component of the 
central spin, a number of non-trivial effects are known to appear, despite its apparent 
simplicity. One instance is the crossover from overdamped to underdamped behaviour 
as the coupling strength increases, similar to the spin-boson model where the TLS is 
coupled to an oscillator bath [H [2] . For a more complete account on the dynamics of a 
spin coupled to a single bath we refer to Refs. [HI [T6] . 



3. Time evolution 



In the absence of bath dynamics the x- and y-components of the total spin of the 
respective environments, Jfo^ = J2k K J^^t = 'Ylk "^fc ' conserved quantities. We 
write the total Hilbert space as a tensor product Tis ® '^^b- The total 2^^ dimensional 
Hilbert space of the baths Tie decomposes into invariant subspaces Tij^^'™^-'. These 
are labeled by the eigenvalues mi, m2 of the total spin operators /f^t and J^ot- Each 
rrij {i = 1,2) runs from —N/2 to N/2. The Hamiltonian ([1]) acts on the subspace 

7^s®H^^""'^ as 

H\a)\ai,a2;mi,m2) = (cuqS"' + S^'migi + S^m2g2)\cr)\ai,a2;mi,m2) ,(4) 

for \ai,a2',mi,m2) G T-^^i'™^^ and |cr) G Tis- Here a,, which is not important in 
the following, labels the irreducible representation. We therefore can write H most 
conveniently as a direct sum H = ®mim.2=-N/2^™-i^2 where 

= S^-fi„,^, ®Ia„,^ ®Ia™^ . (5) 
Here Ia is the A x A unit matrix and the parameters Am^, 

^- = ( 7v/- ) • '=''^ 

measure the dimension of the invariant subspace T-L^^'"^^\ From Eq. ([5D it is clear that 
the effect of the environment is to give rise to an effective magnetic field flmim2 = 
(mif/i, m25'2, i-^o)- However, this effective magnetic field is different from the static 
magnetic field pointing in the 2;-direction, since it does not take a single value but 
rather is a distribution characterized by the degeneracy coefficients Am^. 

Since -ffmun2 acts on the subspace of the environment, we will often write Hmim2 = 
S ■ ^mim2 for short. The eigenvalues of Hmim2 are ±^^17712/2, where we introduced the 
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frequency flmim2 = i^o + S'l'^i + 92^2)^^'^ ■ W^e denote the eigenstates of -ffmim2 by 
|±,mi,m2). They are related to the eigenstates | t), | |) of the non-interacting system 
Hamiltonian Hs by the unitary transformation 

|+,mi,m2) j ^ [ cos0rmm2 sin6'„,^2e 

|-,mi,m2) J y -sin6'„,„2e-^'^'"i™2 cos6'„^„2 

where the angles 4>mrm2 and 6'mim2 are given by 



Km-i mo 




(7) 



4>mim2 = arctan ( ) (8) 

2 n _ ^mima + ^0 /qv 

^mim2 

The case of a single environment is recovered by setting g2 or, equivalently, 4>mim2 equal 
to zero. 

The time evolution operator is straightforwardly derived from Eq. ([5]). We obtain 

U = ®Zlm2=-N/2Urmm2 with 

Umim2it) = COS ('l^^mima) h + ^i (2 ""'"^ g ■ (^^^^^ ^ (10) 

We can decompose an arbitrary system (central spin) operator O as O = 
®m(^m2=-Af/2^™i™2- particular we are interested in the Heisenberg spin operator 
S{t) and its commutators and anticommutators at different times. Using (ITU]) , we find 

Smim2it) = COs{Q„iim2t)S{0) - Sm^Qrmmit) {S (0) X n^^mj 

+ [1 - COS(1]„1„2^)](5(0) ■ n^im2)^mim2 (H) 

-^[^^,^2(^)' 5;^i,n2(0)] = cosin„,,^,t)e,,,S\0) (12) 

+ sin(fi^j^2t) (^5ijn^i„^2'S^(0) - S\0)nl^^„ 
+ [1 - cos(^]™,^2t)]n^^„^ (^S^(O) X n^rm^j , 

2{*S'm^m2 (t), iS^^j^^ (0)} ~ COs(i7m^m2'^)'^ii ~ Sin(fim.im2^)£«jfc'^mim2 

+ [1 - cos(^]^,^2t)] n^^^^n^^^^^ . (13) 
The vector nmim2 is a unit vector pointing in the direction of the effective magnetic field 

^mim2 ■ 

If the density matrix p{t) of the total system is initially invariant under rotations 
within a subspace T^j^^'"'^'' due to the trivial action of the Hamiltonian in this subspace, 
this invariance will persist at all times. In particular, the density matrix can be written 
for all times as p{t) = ®m^m2=-N/2P^i"i2it)- This means that p{t) shares for all times 
the block structure of the Hamiltonian. If p(0) fulfills this condition, the expectation 
value of an arbitrary system operator 0(t) with respect to p(0) can then be written as 

(O(t)) = tr[p(0)O(t)] 

N/2 

^ ^ '^mi -^m2 ^»nim2 (^)Prrtim2 (0) ("^^) 
mi ,m2=—N/2 
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Figure 1. The magnetization of the TLS is plotted with time for various values of the 
system bath interaction, in zero-field. In the presence of two different baths a change 
of sign arises during the time evolution. The sign reversal reaches its maximum value 
when the two baths are identical. The total number of spins in each bath is taken 
to be iV = 100. The system bath couplings given in the inset are dimcnsionless 

i9^/V9^+9i)■ 



In the following we will analyse the expectation values of the operators ffTTj) to f|T3l) 
with respect to an initially unpolarized bath. Since the magnetic field applied along the 
^-direction only affects the central spin, we take this to be initially in the ground state 
determined by Hs and consequently choose the initial density matrix as 

p(^) = ^1 T)(T I®i2-®V • (15) 

We immediately see that, in this state, (5'*(0)) = —5^/2. 

Using the above formalism for evaluating the dynamical properties of the TLS 
we shall now calculate various quantities for the TLS operators and study the effects 
brought about by the coupling to two different baths. 



3.1. Expectation values 

The expectation values {S^{t)) can now be calculated using Eqs. (fTTi) and (fTSjl . We find 

_1 

nil ,1^12=— N/2 

— COS{Qmim2t) + -^^ 1 (16) 

mim2 mim2 / 

{s^{t)) = {sy{t)) = (17) 



Quantum frustration of dissipation by a spin bath 



8 




Figure 2. Same as Fig. [T] but for a nonzero field sucli that loqI \/ g\ g\ = 2. In 
the presence of the external field the magnetization for the single bath case decays 
non-monotonically and saturates to a non-zero value. 



In the absence of external field ujq = 0, we would expect that the initial polarization of 
the TLS would decay faster in comparison to the single bath case, since the total number 
of spins with which the TLS is interacting is doubled. In Fig. [1] we have plotted the 
time variation of the {S^{t)) for various values of gi,g2 keeping g = \/ g\^ g\ constant. 
As one can see, for a single bath the polarization decays to zero very fast, whereas in 
the presence of the second bath, the decay is comparatively slow. In contrast to the 
case of single bath, one observes a change of sign in the time-dependent behaviour of 
the polarization of the TLS, indicating the presence of a nonzero field. 

In Fig. [2]the time variation of the {S^{t)) is plotted for non-vanishing external field 
Uq 7^ 0. We find that the polarization saturates to a finite value at long times with 
faster oscillations in the case of a symmetric double bath, which is consistent with the 
spectral properties discussed later in the text. Inspection of Figs. [1] and [2] reveals that 
the change of sign occurs in the single bath case only if there is a nonzero field, while it 
is observed for both zero and nonzero field in the case of a symmetric double bath. Since 
the baths are completely unpolarized (peaked at m = 0), it is clear that the effective 
field responsible for this change of sign can only stem from the competing effect of two 
baths coupled to non-commuting components of the central spin. 

3.2. Quantum purity 

For an arbitrary density matrix p, purity is defined asV = Trp^. Purity is a convenient, 
basis-independent measure of the degree of coherence, if p is the reduced density matrix 
of the central spin (that which results from tracing out the bath degrees of freedom 
in the total density matrix). In our case the decay of purity is directly related to the 
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Figure 3. The purity of TLS is plotted against time for various values of the system 
bath interaction. Though the initial decay rate strongly depends on the combined 
interaction strength \/ g\+ g\ (which is same for all the three cases), the TLS looses 
its initial polarization faster in the presence of a single bath in comparison to the case 
of two baths. The external field is set to zero and the number of spins in each bath 
are taken to be equal with N = 100. The system bath couplings given in the inset are 
dimensionless {gij \J g\ + g\)- 



relaxation of the spin expectation values to equilibrium. 

m = l + f2(s\t))'. (18) 
1=1 

The result is plotted in Fig. [3l We notice that, both in the single and double 
symmetric bath cases, the purity decays fast to its minimum value 1/2. In the symmetric 
case, we notice a small, short-lived revival that may be interpreted as a weak form of 
decoherence frustration which however does not affect the long time behavior of the 
central spin. 

4. Density of states 

To get further analytical insight we introduce the function 

^ -7V/2 

-P(^) = 22Af+l >^mi>^m2[^{^ - ^mimi) + ^{^ + ^mimi)] (19) 

mi ,m2=—N/2 

which is essentially the density of states, normalised to fulfill the sum rule J duoDiuo) = 1. 
The expectation value (S'^(t)) is related to D{u) by 

{S\t)) = dujDiuj) + • (20) 
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Figure 4. Density of states as a function of frequency in normalized units. 
D{uj) can further be evaluated by using the approximation 



2N ^ - V TlN 



dme-^""'/^ , (21) 



for the binomials A^i, which is known as Laplace-de Moivre formula in probability 
theory [19] and which is valid only for large N. We find 

iV^72sin(2^/'° V iV^2sin(2^.) ) ^\ W^^Wr 
where the total coupling strength g and the angle 9g are defined by 

gi = g cos dg 

g2 = gsmeg, (23) 

and Iq is the modified Bessel function. The case Og = corresponds to the single bath 
and the case 6g = 7r/4 corresponds to two identical baths. In Eq. (122!) the limit of a single 
bath can be taken by using the asymptotic expansion of the modified Bessel function 
lim^^oo Io{z) = /y/2Tcz + . . . The results are plotted in Fig. HI As the coupling of the 
TLS with the baths becomes symmetric, i.e. 6g —>■ it /A, the density of states peaks at 
a frequency away from ujq. As 6^^ — > 0, this peaks shifts towards t^o, which is expected 
for the case of TLS coupling to a single bath. Thus the density of states can by itself 
reveal the frustrating effects of decoherence more elegantly. 

For the particular case of Hs = {ujq = 0), one can obtain simplified expressions 
for the density of states. We obtain 

Dico) ^ ^e-V^'^ , g, = g,g, = 
gVN 

D{u)c^^e--'/^'\ g, = g2 = g (24) 
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In the case of a single bath the peak is at = 0, where as for two baths, the peak is 
shifted to = a/ N/2g. We note in this respect that, if the two baths were coupled to 
same spin component, then the behavior would be similar to that of an effective single- 
bath coupled to one spin component. In such a case the peak in D{uj) would remain 
at = 0. Thus, the emergency of a peak at a; 7^ may be viewed as a frustration of 
dissipation due to the competition between two environments coupled to non-commuting 
spin components. 

We end by noting that D{uj) is the Fourier transform of {S^{t)S^{0)) and is a 
measure of the strength of the transitions induced by a periodic perturbation ~ cos (cut). 



5. Correlation functions 

We now investigate the spin-spin correlation functions defined by Cij{t) = (5'*(t)5'-' (0)). 
Specifically, we focus on its symmetrized and antisymmetrized versions, Sij{t) = 
i({5*(t),5^(0)}) and A,{t) ^ -t{[S^{t),S^{0)]). 

Since the system is initially in an eigenstate of S^, the symmetrized autocorrelation 
function in z-direction is simply S^zit) = —^{Sz(t)). Using the general formulae of 
Sec. [31 we find for the transversal symmetrized auto-correlation functions 

^ ^ / ijj^ — |— ^^^TTl^ q'^TTV^ \ 

'Sxxit) = r,2N+2 -^miAma ( 7^^ COs{Hrnrm2t) + 1 (25) 

mi,m2=-iV/2 ^ "^i™2 ^Vima/ 

and a similar result is obtained for Syy{t) by exchanging indices mi and m2. AH 
symmetrized cross-correlation functions Sij{t), {i ^ j) are zero. 

We look at the autocorrelation function in 2;-direction in the case that there is no 
external magnetic field applied ljq = 0. Since S^zit) is proportional {S^{t)) we can use 
the integral representation fl2Ul) . In general, i. e. for intermediate values of 6g and for 
non-zero frequency ujq, the integral (120|) becomes quite difficult and cannot be solved 
analytically. However in some limits closed expressions can be derived. For c^o = we 
obtain 

Szzit) = expi-Ngh'/S) , = 

5..(t) = l-y^texp(-iV/tV8)Erfi(^^j , 0, = 7r/4 (26) 

The antisymmetrized correlation functions are related to the dynamical 
susceptibilities, defined as 

M^)= I ^e-M.,(t). (27) 

In particular, the imaginary part of the susceptibility can be used as a measure of 
the energy dissipated from the system to the bath. Using Eqs. (fT^ and (fUl) we find 
AUt) = 0, 

-1 UJo 
•Axx{t) = r)2N+l ^ ^ ^rni^m2Q Sin(f2mim2^) 

mi,m2=—N/2 
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dujD(uj) — sm(ujt) (28) 
to 

and Ayy(t) = Axxit), where D^u) has been analyzed in the previous section. All anti- 
symmetrized cross-correlation function but Axy{t) are zero. For Axy{t) we find 

_1 

Axy(t^ 22JV+1 ^ ^ ^mi^m2 *-'Os(r2mim2^) 

mi ,m2=~N/2 

dujD{uj) cos(cut) (29) 

In the second lines of Eqs. f l28|) and fl29l) we used an integral representation in terms of 
D{uj). In this form the dynamical susceptibilities are readily evaluated 

Xxx{uJ) = — duj — — 30 

27r y tu^ — uj'^ + sgn(co')20+ 

Splitting Xiji.^) ill its real and its imaginary part, x«j(^) = x'ijil^) + ^Xij{l^)i obtains 
the relation ([22]) 

xL(^) = ^^M, (31) 

which holds for uj > ujq. Moreover x'xxi'^) = for |ti;| < c^o, and xLl""^) = ~Xxx{'^)- 
We can use the approximation fl2Tl) and obtain 

x4M^|^e.p(-lM). (32) 

From the above equations it can be seen that there is a strong singularity at uj = ujq/ui 
addition to the Gaussian spread arising due to the interaction with the bath. For the 
symmetric coupling, this singularity is removed and only a Gaussian spread peaked at 
UJ = ujq remains. Both functions are peaked aXuo = uoq, and hence one can say that, since 
there is no peak shifting there is no frustration. If we try to remove the singularity for 
the single bath case by multiplying \/ iJ^ — lJ^ with yfl-J^iS) then one can immediately see 
that the peak for xfl^ioS) is shifted away from uj = ujq for the symmetric case. Scalings of 
such kind can be avoided by considering other kinds of distributions for the bath spins. 
In the next section we consider bath spin distributions Am with a Gaussian cutoff. 



6. Tailoring the density of states 

In Refs. [3111] respectively [5], [8] similar expressions were obtained for yfl-JyoS) in the first 
case and for the Fourier transform of the antisymmetrized position-position correlation 
function in the second case. In both cases the function under consideration is essentially 
a Lorentzian 
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where Uq is the renormahzed frequency of the system (Larmor frequency, respectively 
oscillator frequency), g^s is the effective damping coefficient. For the detailed expressions 
of uJo, Qes and Z see Eq. (26) of Ref. |lj and Eq. (20) of Ref. The Lorentzian form, 
and in particular the linear behavior for small values of u is typical for Ohmic type of 
dissipation. 

We model our system to best mimic Ohmic behavior. Since Xxxi^) is related to 
the density of states D{uj) by the simple relation (|3T|) . we can focus directly on D{uj) as 
defined in Eq. (fT9|) . We now use a general Am of the form 

A(") - 1 (1) exp (-^J!^) (35) 
"2^^ r((l + a)/2)^l, n)- 

For small values of m the behavior of is dominated by the power m° with a 
characteristic exponent a. We note that the case a = 1 cannot directly be identified 
with an Ohmic bath. In Eq. fl35|) we took a Gaussian cutoff for large values of uj with 



a cutoff frequency chosen as a/ N/2 in order to make contact with the former results. 
It is a well known fact in the theory of open quantum systems that the specific form 
of the cutoff function is not relevant [2]. The function A^, is normalized such that 
/ dmX{m) = 1 holds. We see that the form of A^ described in Eq. (jS]) is just a special 
case of Eq. fl35|) corresponding to a = [see also Eq. fl2Tl) ]. 

In a calculation which is similar to that performed in Sec. 4, we obtain for the 
density of states 

= r((a+ lV2Wsm(2(),) (iV9^cos(29,)) ' (-^ -or'^«(— o) 

2 , ,2\ \ / /I A ,2 , ,2\ 



/ 4cot2gg(a;^-^g) \ / Ajuj - . 
\ Ng^ sm{2e,) ) "''P \ Ng^ sin' (29,) ) ' ^^^^ 

where, as before, g'^ = gf + g2 and the angle 9g is defined in Eq. f l23|) . Moreover we 
have introduced the modified Bessel function 1^/2 of order a/2. For the transverse 
susceptibility we find in the two limiting cases Og = and 9g = 7r/4 the expressions 

1/2 2\Q^ / A \ Q+1 



r(a + i) 

1 iu,' - ^'f-'>" ^0 





71 

4 


2 ' 




Ng\ 


^ (a-l)/2 


0,= 


. 



^xxv ; 2Ar^2 r((a + l)/2) 

which satisfies the general properties given after Eq. (1311) . The transverse susceptibility 
is zero for uj = uq. This zero value is found because the distributions A^ is centred at 
m = 0. If the distribution is shifted to be centred at a non-zero m value, then the value 
of x"(a;)i2 will be non-zero a.t u = uq for sufficiently large a. In Fig. [5] we have plotted 
the transverse susceptibility for the symmetric {9g = 71/4) and single-bath {9g = 0) 
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Figure 5. Plot of the transverse suceptibility Xxxi'^) for an Ohmic type of density 
of states a — 1 and for three different angles 6g = 0,7r/8, 7r/4. The other parameters 
are luq = 10, g = 1 and N — 1000. The offset at which the function is normalized is 
e = 0.1 



cases, as well as for an intermediate situation. In order to compare with other results, 
in particular with the curves obtained in Ref. we have normalized Xxxi^) so that 
Xxxi^o + e) = 1) where e is a convenient small offset which is chosen for proper scaling 
and comparison. We note that making e too small shoots the peak to infinity in those 
cases where lim^^o Xxxl'^o + e) = 0. Though the natural sum rule / dco uJXxxi^) = '^o 
is spoiled by adding this epsilon the essential physics behind is unaffected. When there 
is no interaction with the bath, x"{^) is a delta function peaked at ci; = uq. In the 
presence of one bath the peak broadens with the maximum still located at u = Uq. 

Surprisingly, the peak at = cuq disappears when we introduce a second bath which 
couples to a different component of the central spin. A similar shift was reported in 
Refs. [31 H]. It results from a pure frustration effect due to the non-commuting nature 
of the spin operators. 

In Fig. [6] we have plotted the transverse susceptibility for three different types 
of infrared behavior a = 0,1, 2, mimicking a subohmic, an Ohmic, respectively a 
superohmic bath. Dissipation decreases as the power a increases, as expected for general 
dissipative quantum systems p. On the other hand we see that the frustration effect 
of an additional bath increases with increasing power a (not shown) . 

7. Conclusions 

We have analyzed the spectral properties of a two-level-system which is coupled to one 
or two dissipative baths through non-commuting observables. A peak in the spectrum 
at a nonzero frequency reveals the existence of an effective magnetic field experienced by 
the central spin. We have seen seen that the coupling to a second bath enhances rather 
than diminishes that effective field and, with it, the dynamics of the central spin. In 
the extreme case of a zero external field, the nonzero field is generated solely from the 
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Figure 6. Plot of the transverse suceptibility x'xxi'^) f'^'' the symmetric case 9g ~ it/ A 
for the values a = (dotted), a = \ (dashed) and a — 2 (full). The parameters are 
cjQ = 10, g = 1 and N = 1000. The offset at which the function is normalized is e = 0.5 

competition of two environments coupling to non-commuting spin components. This 
fact is remarkable if one notes that the baths are assumed to be initially unpolarized. 

These physical effects arising from the non-commuting nature of the spin operators 
are a general feature of the dynamics which does not depend on details such as the 
Markovian or non-Markovian character of the reduced system dynamics or the strength 
of the system-bath interactions. In fact, we explicitly prove that the emergence of a 
peak in the spectral function is compatible with a fast decay of the quantum impurity. 
This suggests that, while dissipation is inhibited by the competition between the two 
baths, decoherence is not, at least in the long time behavior. Our present results are 
entirely consistent with the results of Refs. [71 [H] for a harmonic oscillator (equivalent 
to a large or quasiclassical spin), where a fast decay of the quantum purity was found 
to coexist with a weak form of suppression of dissipation. Here we have proved that the 
phenomenon of frustration of dissipation also exists for a dissipative two-level-system 
which, given its reduced dimensionality, is much more quantum in nature than the 
harmonic oscillator. 

The intuitive idea that two competing environments attempt to measure non- 
commuting observables and thus fail to generate decoherence sounds appealing but 
may be misleading. The statement would be true if the only possible result of 
a quantum measurement were to select a narrow distribution of eigenstates of the 
measured observable, since two non-commuting observables cannot be simultaneously 
well defined. However, a possible outcome of the coupling to a dissipative environment 
is that the reduced density matrix, while becoming diagonal in the representation of the 
eigenstates of the measured observable, may display a broad probability distribution in 
that representation. In the limit in which that distribution is very broad, the reduced 
density approaches the identity matrix, which is invariant under a change of basis. Thus 
a reduced density matrix may be simultaneously diagonal in the representations of two 
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non-commuting observables, provided it is close to the identity matrix. This is what 
actually happens to our central spin-i, as is clearly revealed by the quantum purity 
tending to its minimum value 1/2 at long times, both for a single and a symmetric 
double bath. As we have seen, this feature is compatible with the reinforcement of the 
central spin dynamics resulting from the competition of the two environments. The 
upshot of the present study on the effect of competing environments is that, at least for 
dissipative two-level-systems, it may be misleading to speak of quantum frustration of 
decoherence and is more appropriate to introduce the concept of quantum frustration 
of dissipation. 

We cannot rule out however the possibility that the coexistence of decoherence and 
dynamics enhancement by two competing environments is a particular property of our 
dissipative model where the bath has no internal dynamics. A firmer conclusion on 
the existence or absence of decoherence frustration will require an understanding of the 
behavior of genuinely quantum properties such as purity or pair entanglement in the 
presence of competing environments with internal dynamics. 
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